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ABSTRACT 

We study temperature and density sensitivities of ratios of Si XI soft X-ray 
emission lines, in the wavelength range of 43 — 54A . The typical temperature of 
the formation of the analyzed lines is around 1.6 MK, which makes this analysis 
complementary to the analysis of He-like triplets being sensitive to hotter plasma. 
We present theoretical calculations and compare them with ratios obtained from 
high-resolution X-ray spectra of five solar-like stars: Procyon, a Cen A&B, e Eri, 
and Capella. We find that our results are in good agreement with results obtained 
by other authors through different diagnostics, namely the analysis of density- 
and temperature-sensitive He-like triplet lines. We further estimate the coronal 
pressure and filling factors from Si XI lines in this study. 

Subject headings: atomic data -techniques: spectroscopic - stars : coronae - 
stars: late-type - X-rays : stars 

1. Introduction 

Determinations of elemental abundances and structures in a various of stars belong to 
the most important tasks of observational astronomy and the understanding of the evolu- 
tion of chemical elements. The temperature structure of emitting layers such as the stellar 
coronae must be known before elemental abundances can be determined. The density is 
another important physical parameter to describe the magnetically confined plasma in outer 
atmospheres of late-type stars, which can be used to estimate the spatial information for the 
stellar coronae. For those pre-main sequence stars, the constrained density can give us clues 
for the origin of the X-ray production such as magnetics alike the stellar coronae or a shock 
due to accretion processes. 



- 2 - 



Before the launch of missions with high resolution spectrograph such as the Chandra 
and XMM- Newton, a determination of the electron temperature of X-ray emitting layers 
usually adopts a technique of global-fitting to spectra. The determined temperature struc- 
ture from this technique strongly depends on the determination of continuum level of the 
observed spectra. The measurement of electron density is not possible for non-solar stars 
with previous X-ray satellites such as ASCA and ROSAT configured with low- and medium- 
resolution spectrographs. The observation data provided by high-resolution spectrograph 
on board new generation satellites Chandra and XMM- Newton, allow us to estimate this 
parameter. Individual emission lines can be resolved from high-quality spectra, further the 
electron temperature can be derived from line intensity ratios such as G value (between inter- 
combination % and forbidden / lines vs resonance line r) of He-like ions which was described 
by Porquet et al. (2001) in detail. Since then, the estimation of the electron density becomes 
possible from X-ray emission line ratios such as R value between % and / lines He-like ions. 
Using observation data with lower spectral resolution, n e and emission measures EM can 
not be obtained independently, such that no emitting volumes V could be estimated from 
the relation of EM = n\ V. The structural information of the outmost atmosphere of stars 
was accessible from the analysis of lightcurves for a few special systems such as rotational 
modulation, eclipse mapping, etc. However, the solution of the reconstructed differential 
EM is always nonunique, and involves the difficulty of disentangling time-dependent versus 
geometric variability (Schmitt et al. 1996; Schmitt & Favata 1999; Giidel 2004). Since the 
operation of the two satellites with high-resolution and high-collecting areas, the two param- 
eters have been derived from spectroscopic methods for every class of X-ray sources by many 
observers (Brinkman et al. 2000; Canizares et al. 2000; Kastner et al. 2002; Ness et al. 2002, 
2003; Raassen et al. 2003; Stelzer & Schmitt 2004; Schmitt et al. 2005, etc.). Using He-like 
triplet lines, a systematic investigation of the electron temperature and density for late-type 
stars with various levels of activity, has been made by Ness et al. (2002, 2004) and Testa et 
al. (2004). However, the UV field significantly affects the He-like triplet lines. Recently, Ness 
& Schmitt (2005) presented a density diagnostic ratio being independent from the UV field, 
that therefore allows to disentangle n e and UV field that both affect the O VII triplet lines. 
Additionally, the derived density has also been used to estimate the production mechanism 
(alike stellar corona or accretion) of the X-ray emission in pre-main sequence stars (Stelzer 
& Schmitt 2004; Schmitt et al. 2005). 

Soft X-ray emissions of L-shell ions of silicon have been detected extensively in spectra 
of stellar coronae. These emissions carried out information of the two physical parameters 
for the solar-like coronae. In our previous work, we investigated the emission lines of Si X, 
and diagnosed the coronal electron density using them for late-type star — Procyon (Liang et 
al. 2006a,b). The soft X-ray emissions arising from 3d — 2p, 3s — 2p and 3p — 2s transitions 
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of Si XI have also been detected obviously for some cool stars, such as Procyon and a Cen 
A&B (Raassen et al. 2002, 2003). And the temperature (1.6 MK) of maximum formation 
of Si XI ion in the ionization equilibrium (Mazzotta et al. 1998), is the typical temperature 
of corona of stars with low activities, which makes the analysis of Si XI soft X-ray emission 
lines interesting. 

EUV features arising from 2s 2 -2s2p and 2s 2 -2p 2 transitions of Si XI have been used to 
derive the electron temperature and density for solar flares. Keenan et al. (1995) presented 
seven line ratios whereas these ratios are usually sensitive to both the electron temperature 
and density of the emitting plasma. Moreover measured ratios obtained with the Naval 
Research Laboratory's S082A spectrograph on board Skylab, are larger than theoretical 
high-temperature and high-density limits (Keenan et al. 1995). Later, Lang et al. (2001) 
calculated these line ratios again and compared with solar observations obtained with normal 
incidence spectrometer on board SOHO CDS, as well as laboratory measurements performed 
by Konig et al. (1996). In contrast to the EUV lines, there has been little work on the soft 
X-ray transitions of Si XI in astrophysical literatures, probably as a result of the limited 
availability of high quality spectra of stellar in this wavelength region. 

This paper is organized as following, we firstly investigate the Si XI soft X-ray emis- 
sion lines and present theoretical calculations of temperature- and/or density-sensitive line 
intensity ratios in Sect. 2. A brief description of observations and extractions of line fluxes 
is presented in Sect. 3. Results and discussions are outlined in Sect. 4. Our conclusions are 
given in Sect. 5. 

2. Theoretical line intensity ratios 

Our model of Si XI consists of energetically lowest 350 fine-structure energy levels 
belonging to configurations 2s 2 , 2s2p, 2p 2 , 2131', 2141', 2151', 2161' and 2171' (l=s, p; 1' = 
0, 1, ...,n — 1). Some energy levels are replaced by available experimental values from the 
NIST database 1 and results of Coutinho & Trigueiros (2001). The theoretical calculation is 
performed using a fully relativistic method — flexible atomic code (FAC) developed by Gu 
(2003, 2004). Self-consistent results of electron impact excitation, de-excitation rates and 
Einstein ^-coefficient (including electric- and magnetic-dipole, as well as quadrupole transi- 
tions) among the 350 levels were also included in present model. The input atomic data of 
our model significantly differs from that of CHIANTI code, in which only 46 low-lying energy 
levels, and the excitation and radiative decay among the 46 levels were included. Present 



1 http://physics.nist.gov/cgi-bin/AtData/main_asd 
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model is an extensive one, and cascade effects on the upper levels of interested transition 
lines from higher excited levels up to n = 7 (n is the main quantum number), have been 
considered. In our recent work (Liang et al. 2006c), a detailed assessment for the calculated 
data is performed, which reveals that the atomic data is reliable for diagnostic application. 
Here, a brief comparison of our data with the input atomic data of CHIANTI code is made 
as shown in Fig. 1. The top panel of Fig. 1 illustrates the comparison of different theoret- 
ical predictions of energy levels with the available experimental ones (the NIST data and 
that of Coutinho & Trigueiros (2001)). This panel indicates that our energy levels agree 
with the experimental values within 0.7%, which is better than the theoretical data in the 
CHIANTI database. NIST only lists the energy levels being less 20 Ryd. For higher energy 
levels of above ~20Ryd, Coutinho & Trigueiros (2001) reported some values and determined 
the energy levels from observed wavelengths by an interactive optimization procedure using 
the program ELCALC (Radziemski & Kaufman 1969). Comparing with their energy levels 
also can benchmark our prediction in a certain extent. A better agreement is shown by the 
filled circles in the upper panel of Fig. 1. The comparison of weighted oscillator strengths 
is illustrated in the middle panel of Fig. 1, which reveals the two different data is com- 
parable within 20% for most strong transitions (#/>0.1). For those transitions involved in 
below interested ratios (filled up-triangle in Fig. 1), the two different predictions show a 
good agreement within 20%. The prominent 3d-2p transitions exhibit a better agreement. 
For some weak transitions (gf <0.1), large discrepancies up to several factors appear. We 
suggest such differences are due to different inclusions of configuration interaction (CI) in 
the different data sources. The work of Aggarwal et al. (2005) and Zeng et al. (2005) has 
revealed that the consideration of an elaborate CI is very necessary, which support our per- 
formance of atomic data calculation. For the effective collision strength (T) obtained under 
the Maxwellian distribution, only the data at a temperature of 1.0 MK (the bottom panel of 
Fig. 1) is compared for concision, because the temperature is the typical coronal temperature 
of stars with less activity. 

Proton impact excitation is important process for population and de-population of the 
lowest levels with n = 2 complexes. Data from Ryans et al. (1998) was employed in present 
model, in which the close-coupled impact-parameter method was adopted in the calculation. 

The theoretical line intensities are obtained by solving steady-state rate equations with 
optically thin assumptions, as described by Gu (2003) and Behar et al. (2001). Such assump- 
tions are based on the following reasons. Firstly, photon-excitation and de-excitation rates 
are negligible in comparison with the electron impact excitation rates in stellar corona, toka- 
mak and laser-produced plasmas, so the processes aren't considered; secondly, the rates of 
ionization to and recombination from other ionic stages are lower than bound-bound transi- 
tion rates; finally, emitted photon directly escapes the plasma, does not collide with particles 



- 5 - 



again. The line intensities are calculated over a wide range of the electron temperature and 
density. Here, we pay careful attention on five prominent lines of Si XI, as listed in Table 1. 
The 3s -2p line at 52.306A is a strong and isolated line in these selected features. A detailed 
analysis indicates that following emission line intensity ratios of Si XI ion are sensitive to 
the electron temperature and/or density: 



Ri 


= /(52.306 A)//(43.743 A), 


R2 


= /(52.306 A)//(49.207 A), 


R 3 


= 7(49.207 A)/7(46.39lA), 




= 1(52.306 A)//(46.391 A), 


R 5 


= 7(43.743 A)/7(46.391 A), 


R% 


= 7(46.391 A)/7(46.283 A), 



and 

R 7 = 7(52.306 A)/7(46.283 A). 

Therefore they may be useful diagnostic techniques for the electron temperature and density 
for all kinds of hot plasma. 

Fig. 2 plots these ratios as a function of the electron temperature at n c =1.0xl0 9 cm -3 , 
and vs the electron density at logT c (K)=6.2 which is the temperature of maximum formation 
of Si XI(see solid line in the figure). The ratios are given in photon units, and used throughout 
the paper. The ratio R± varies from 1.7 to 0.9 over the temperature range of logT c (K)=6.0— 
6.4, and no variation occurs through a large range of density 10 7 — 10 12 cm~ 3 , which is the 
typical condition of astrophysical, tokamak and EBIT plasmas. This indicates the ratio Ri 
is a powerful diagnostic technique of the electron temperature for the hot astrophysical and 
laboratory plasmas. Ratio 7? 2 exhibits the same behavior as Ri, yet this ratio shows a less 
temperature-sensitivity than that of R±. 

The ratio 7? 3 appears to be sensitive to the electron density (shown in Fig. 3), simulta- 
neously it displays a certain sensitivity to the electron temperature. For example, increasing 
logT c (K) from 6.2 to 6.4 leads to 8 per cent variation at n c =1.0xl0 9 cm" 3 . The feature at 
46.391 A being correlative to ratios i?3, R4 and R§ has been identified to be 46.410 A line 
arising from 2s2p 3 P2 — > 2s3d 3 D3 transition of Si XI in work of Raassen et al. (2002). After 
adjusting by experimental value from the NIST database, a better wavelength 46.401 A is 
obtained in present work. We noticed that the emissivity of the line at 46.404 A is about 
10 per cent relative to that of 46.401 A, whereas the two lines can not be resolved using 
present available spectrographs. The expected blending of 46.404 A line is rather smaller 
than both the theoretical uncertainty and statistical error in the line flux measurement, 
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however the blending effect was considered in this study. Ratio R 4 is also sensitive to the 
electron density, but it is insensitive to the electron temperature in the density-sensitive 
region, as shown in the lower panel of Fig. 2. The variation is less than 4 per cent when 
logT c (K) changes 0.2 dex. This indicates R4 is a powerful n e -diagnostic tool for the line 
formation region. As R%, the blending effect in 46.391 A line has been also considered in the 
prediction of R4. Ratios R$ — R7 show obvious n c and T c -dependence, as shown by curves of 
i? 5 and R 6 in Fig. 2. For clearness, the curve of R? vs T c is not given, only its dependence 
on n c is plotted in Fig. 3. The resolved lines and their blending contributions are listed 
in Table 1. From above analysis, we distinguish that R± and i? 4 are powerful T e and n c 
diagnostic method, respectively. Using the properties of Si XI, we can estimate the electron 
temperature and density in Si XI emitting region for cool stars. 



3. Observations and line fluxes 

In this work, we re-analyze the spectra of stars including three normal dwarf stars, 
i.e., a Cen A&B, Procyon and e Eri, and an active binary system Capella. The detailed 
descriptions and analyses for these stars are available from many literatures (Schmitt et al. 
1996; Mewe et al. 2001; Audard et al. 2001; Ness et al. 2001; Raassen et al. 2002, 2003, etc.). 
Here, we briefly describe the observations and data reductions for these stars, a Cen A 
and B are firstly resolved from spectral observation with Low Energy Transmission Grating 
Spectrometer (LETGS) instrument on board Chandra observatory, a detailed description 
about it was given by Raassen et al. (2003). All observations use the LETGS combined with 
High Resolution Camera (HRC) detector. Reduction of the LETGS datasets adopts CIA03.2 
software with the science threads for LETGS/HRC-S observations. For the extraction of the 
spectra of a Cen A&B, a similar procedure used by Raassen et al. (2003) is adopted here 
under the CIAO environment. 

Line fluxes are obtained by locally fitting the emission features using multiple overlap- 
ping Gaussian profiles, together with a constant value representing background and (pseudo- 
continuum emission, which was determined in line-free regions such as (51 — 52A). The pos- 
itive order spectra was used to derive the line fluxes because some lines such as 52.306 A of 
Si XI lie in the gap between chips in the negative spectra. Determinations of the line fluxes 
of Si XI lines are difficult because a large amount of emissions of highly ionized magnesium, 
sulfur and argon is also present in this wavelength range of 43 — 53A . Here we choose Procyon 
as an example to illustrate. At the shorter wavelength wing of Si XI line at 43.743 A , a weak 
line of S XII at 43.649 A appears as illustrated in the left panel of Fig. 4. Two emission 
lines of S IX at 49.107 A and 49.324 A broaden the shoulder of Si XI line at 49.207 A as 
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shown in the right panel of Fig. 4. By the best-fitting to the spectra, observed line fluxes are 
obtained as listed in Table 1 with la statistical error. In this table, labels such as 3a and 36 
in first column denote this line is blended. For Procyon and a Cen A&B, all prominent lines 
of Si XI have been identified by Raassen et al. (2002, 2003), and line fluxes are consistent 
with those derived by us within la statistical error. However, for e Eri and Capella, only 
the line at 49.207 A is reported in available literatures. 

4. Results and Discussions 

4.1. Determinations of T c 

Sect. 2 has revealed that Ri and R 2 are good T e diagnostic tools. By comparing observed 
line ratios R\ and R 2 with the theoretical predictions as shown in Fig. 5, the electron 
temperature of the Si XI emitting region is obtained, as listed in Table 2. The observed 
ratio Ri of a Cen A slightly higher than other observed ratios, and the value is up to 1.72. 
This implies that the diagnosed temperature is lower than that of other stars. In other cases, 
the observed ratios of R\ are ranging from 1.2 to 1.4. In consequence, the diagnosed electron 
temperatures lie within ~1.5-2.4 MK. These results are consistent with the temperature 
(1.6 MK) of maximum formation of Si XI in the ionization equilibrium (Mazzotta et al. 
1998). The results also agree with the peak temperature of differential emission measure 
(DEM) constructed from a global spectral fitting for normal stars including Procyon and 
a Cen (Raassen et al. 2002, 2003), and low-temperature component of DEM for active 
binary system — Capella (Audard et al. 2001). For a Cen binary system, our results further 
confirms Raassen's conclusion derived from He-like triplet that a Cen B (K1V) component 
is hotter than the a Cen A (G2V) component (Raassen et al. 2003). 

As stated in Sect. 2, R 2 shows a less sensitivity than that of R 1: so measured ratios are 
compatible with the upper limit of the theoretical prediction within la statistical error. In 
the case of Procyon, Raassen et al. (2002) have reported that the Si XI line at 49.207A is 
contaminated by Ar IX line at 49.181 A which can not be resolved out by present LETGS 
instrument on board Chandra observatory. Therefore, the correct extraction of the fraction 
of Ar IX line is very important. Fortunately, Lepson et al. (2003) also detected this line in 
their recent experiment conducted in EBIT measurement for L-shell spectra of highly charged 
argon. Using the experimental ratio (1.11) between the lines at 48.737 A and 49.181 A of Ar 
IX measured in the EBIT facility, we extracted the contribution of Ar IX line and rederived 
the ratio R 2 . The weak ^-dependence of the ratio of the two Ar IX lines supports the direct 
application of the experimental value. Taking Procyon star as an example, the Ar IX line 
at 48.730 A has been obviously detected and no other contamination. We conclude that the 
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Ar IX line at 49.181 A occupies about 18 per cent to the total flux around 49.207A. The 
intensity of the Ar IX line is nearly 2a errors. By considering the ~18 per cent contribution, 
the ratio i?2 increases from 0.52±0.13 to 0.64±0.14, accordingly the derived temperature 
drops from 7.07 to 2.24 MK which shows consistency with that derived from i? 4 within error. 
For our other sample, we measured the line flux of Ar IX around 48.730A, and performed 
similar procedure as in Procyon. The derived temperatures show an good agreement with 
those from Ri within la error. We also notice that the intensity of Ar IX line relative to 
Si XI lines changes with the coronal temperature, and specifically, a lower contamination 
from Ar IX lines is expected for more active stars. 



4.2. Determinations of n e 

Sect. 2 indicated that _R 4 is a powerful diagnostic tool for the electron density in the 
Si XI emitting region of hot plasma. By comparing the observed ratios with theoretical 
prediction, we derived the electron density in the Si XI emitting region of coronal plasma, 
as shown in Fig. 5. Here, the comparison of i? 4 together with R 3 is illustrated in Fig. 6. 
Table 3 lists the observed line ratios with la statistical error and the diagnosed density for 
our sample. 

The first impression of Fig. 6 is the uncertainty of observed line ratios appearing to be 
large, and almost all the observed ratios of i? 4 agreeing with high-density limit in statistical 
error with exception of a Cen B case being agreement with low-density limit. A long ex- 
posure times are necessary to better constrain the density from _R 4 . Yet present line ratio 
measurements still can give us some clues for the Si XI emitting region, that is the lower or 
upper density limit can be constrained from our measurements. We note that mean value 
of the electron density for our sample is above 2xl0 8 cm -3 , which is the typical density of 
solar quiescent coronal plasma. For the ratio _R 3 , the contaminated Si XI line at 49.207A 
is concerted, so the blending effect due to the Ar IX line is considered as above subsection. 
Taking Procyon as an example, R 3 drops from 3.60±1.07 to 2.95±0.86 by taking into ac- 
count the contamination. A density of < 8.8 x 10 10 cm -3 is obtained, which agrees with 
that constrained by _R 4 within la statistical error. Since the large statistical errors, only the 
lower or upper limit of the density is given for our sample as listed in Table 3. 
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4.3. Comparison between present and published results 

The temperature of maximum formation of Si XI is close to that of N VII emission lines, 
therefore the physical conditions in Si XI and N VII emitting region can be compatible. For 
N VII, the temperature can be derived in principle from ratio between Lja line and Lj/3 line. 
However the Ly/3 line generally is faint. So the temperature derived from ratio between Lja 
line and resonance line of He-like ion can trace the this information. Firstly, we compare our 
results of T e with those derived from emissions of H- and He-like N ions as shown in Table 
4. In this table, Lya/r column represents temperature derived from ratio between Lja and 
resonance line r of He-like N VI, while the next column represents temperature constrained 
by ratio between Lyo and sum of triplet lines of He-like N VI, whereas the last column G 
denotes the temperature obtained from the ratio % + f vs resonance line. All the earlier 
results are from work of Ness et al. (2002). Present results are slightly higher than those 
determined from emissions of H- and He-like N, whereas they agree within la statistical 
error. The results from G ratio seem to systematically underestimate T e as reported by 
Testa et al. (2004). 

Secondly, we compare present results of n e with those available from literatures as shown 
in Table 5. In available literatures, the density of cooler plasma of stellar coronae are usually 
derived from the ratio between the intercombination (i) and forbidden (/) lines of He-like 
carbon, nitrogen and oxygen ions. Above comparison has denoted the temperature of Si XI 
(or N VII) emitting region is consistent with that of N VI and N VII mixing region. So 
we presumed the emissions of N VI- VII and Si XI are emitted by the same plasma, that is 
the density derived from different techniques should be comparable. Table 5 reveals present 
results agree with the previous work of Ness et al. (2002) from emissions of H- and He-like 
nitrogen. In the case of a Cen B, the constrained upper limit of density is obviously lower 
than the density determined by Ness et al. (2002) within la error, however they agree within 
2a error, such deviation can attribute to the low S/N. In the work of Ness et al. (2002), the 
effect of radiative field has been considered. 



4.4. Electron pressure and filling factor 

Large errors in ratios of Si XI result in large uncertainties in deduced electron temper- 
atures and densities, which imply large uncertainties in derived electron pressures. Here, 
we adopt the n c corresponding to the best-fit value of R4 to trace the information about 
the pressures and structures in line forming regions. Evaluating now the coronal pressure 
(p = k B n c T c ) for stars studied here as shown in Table 6. For Procyon, present result is 
similar to that derived from N VI (4.4 dyn/cm 2 ) by Ness et al. (2001). Argiroffi et al. (2003) 
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also evaluated the pressure using the derived T e and n e from O VII for Capella, and a up- 
per limit of ~6.0 dyn/cm 2 was obtained, which is in good agreement with present result of 
5.57 dyn/cm 2 . The electron pressure of e Eri appears higher than other cases by a factor 
of ~2, and which is also slightly higher than that (6.60 dyn/cm 2 ) derived by Jordan et al. 
(2001) using UV emission lines of Si III and O IV. 

In terms of these deduced n e and T e , we further estimate the coronal structure through 
filling factors in the line forming region. We firstly derive emitting (coronal) volumes V cor 
according to a formula EM [on = 0.85n 2 V cor , whereas EM ion represents a ion-specific emission 
measure which is determined by the correlation between ion-specific luminosity (Lx,i on ) and 
line peak emissivity (e\(T m )), 

EM lon = %^ . 

Here, the cleanest line at 52.306 A is used to deduce the ion-specific emission measure. 
For the peak emissivity, two different values from present calculation and the Astrophysical 
Plasma Emission Code (APEC) database are adopted. Secondly, we estimate available 
volumes V ava i\ which can potentially be filled with coronal plasma. And an assumption 
is made to derive K,vaii, that is the plasma is confined in an uniform distribution of loop 
structures obeying the loop scaling law of Rosner et al. (1978, hereafter RTV): 



with the electron density n 6i hot and the plasma temperature T hot at the loop top. Giidel 
et al. (1997) regarded the loop-top temperature as equivalent to the hotter component of 
a two-temperature distribution, while the temperature of the hotter component has a such 
relation: 



The density at the loop-top n e ,hot we derive it from the measured density n e (Si XI) assuming 
pressure equilibrium in a given loop, i.e., n e T ion = n^hot^hot- From these considerations, the 
available volumes V ava ii can be expressed as 



1.3 x 10 6 T h 2 



hot 





in cgs units. The ratio of these volumes is defined as the filling factor 



/ = 




- 11 - 



In Table 6, we derive the filling factor for Si XI emitting regions along with the deduced 
electron pressures. A few percent of the filling factor for stellar surface is derived again 
from another line ratios. Present peak emissivity (4.31xl0 _16 phot.cm 3 s~ 1 ) is lower than 
the APEC result (5.72x 10 _16 phot.cm 3 s _1 ) by 25%, which straightforward results in the 
filling factor slightly higher as indicated in Table 6. In our prediction, silicon abundance 
of the solar photospheric of Giidel et al. (1997) is adopted, and the ionization equilibrium 
calculation of Mazzotta et al. (1998) is used to deduce the emissivity for Si XI emission 
lines. The different peak emissivity might be due to the different atomic data performed 
with the different inclusion of CI. Moreover, cascade effects from higher excited levels has 
been considered in present prediction. 

Ness et al. (2004) also derived the filling factor for a large sample of stellar coronae 
using O VII and Ne IX lines, and very small filling factor was found again. Assuming 
balance between radiation losses and the net conductive flux, Sim & Jordan (2003) derived 
area filling factor ranging from the transition to the inner corona for e Eri. The derived 
temperature from Si XI emission is ~2 MK. So the derived filling factors may be comparable 
with ones from O VII. The comparison illustrated in Table 6 confirms this point. Testa et al. 
(2004) also investigated filling factors using Mg XI and O VII for 22 active stars. And they 
concluded that the surface coronal filling factors are ~ 10~ 4 — 10 _1 and ~ 10~ 3 — 1 for the 
hotter plasma with T e ~ 10 7 K and the cooler plasma with T c ~ (2 — 3) x 10 6 K, respectively. 
Using the derived electron temperature and density from Si XI, filling factor of 10~ 3 to 10~ 2 
is estimated from Si XI, which shows an agreement with that derived from O VII by Testa 
et al. (2004) for active stars. 

Since the lower or upper limits of the density are obtained from Si XI for our sample, 
we further present the real constraints on the pressure and filling factors as shown in Table 
7. We estimated the pressure in magnetic loops is larger than ~0.5 dyn/cm 2 for our sample 
except for Procyon. This result is complementary to the work of Argiroffi et al. (2003), who 
derived the upper limits from O VII triplets. The filling factors being less than ~4% are 
estimated, which shows good consistency with results constrained by O VII triplets. For 
a Cen B, the lower limit of / (0.36%) and upper limit of p (2.75 dyn/cm 2 ) are estimated. 

5. Conclusion 

In summary, we investigate L-shell soft X-ray emission lines of Si XI by using a more 
complete kinetic model. In this model, cascade effect from higher excited levels with n = 7 
on the upper levels of the present interested lines has been included. Our analysis reveals 
that some line intensity ratios such as R\ and R 2 , provide useful temperature diagnostics in 
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a temperature range of 5.5 — 6.5 (in logT c (K)), while i?4 is useful to diagnose the density in 
a range of 10 8 ~ 10 cm -3 . 

By comparing the observed ratios with theoretical prediction, the electron temperature 
and density for our sample including normal stars (Procyon and a Cen A&B) and active 
binary system — Capella, have been determined. The constrained temperatures by R\ with 
1.6 — 2.4 MK, are consistent with the temperature of maximum Si XI fractional abundance 
in the ionization equilibrium (Mazzotta et al. 1998) and the peak temperature of the DEM 
constructed from a global fitting of spectra for inactive stars. Since the blending effect in 
i?2 from Ar IX line at 49.181 A , the correct line flux measurement of 49.207 A may induce 
a certain of systemical errors. Using the experimental ratio (1.11) of Ar IX lines between 
48.730 A and 49.181 A, we re-derived R 2 and the temperatures, yet they show a good 
agreement with those constrained by R±. Taking into account the contamination of Ar IX 
line, densities constrained by R3 also show agreement with those by R4 within la statistical 
error. 

Temperatures determined by Si XI emissions slightly higher than ones constrained by 
emissions of H- and He-like N, whereas they are in agreement within la error. So we 
presume the emissions of N VI- VII and Si XI are emitted by the same plasma. The yielded 
temperatures from G ratio are systematically lower than present ones as reported by Testa 
et al. (2004). This may be due to the slightly lower peak temperature (1.4 MK) of the 
contribution function of the N VI triplets. An agreement of the densities from the two 
different line ratios of highly charged silicon and nitrogen, is found. This further suggests 
that the Si XI and N VI- VII emitting plasma may be confined in the same loop. 

A filling factor of 10 -3 to 10~ 2 is derived for our sample from Si XI, and a upper limit of 
4% is obtained. Such results are consistent with the findings of Testa et al. (2004) and Ness 
et al. (2004) that filling factors are much smaller than unity; they also support the findings 
of Testa et al. (2004) that the filling factor of cooler plasma is directly proportional to the 
X-ray surface flux of stars. 

G. Y. thanks Dr. J.-U. Ness, Arizona State University, for his very useful discussion. 
This work was supported by the National Natural Science Foundation under Grant No. 
10433010, 10403007 and 10521001. 
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Table 1. Line fluxes of interested lines for 5 late- type stars, and its identification. The 
theoretical wavelength is from calculation of FAC (Liang et al. 2006c), which is corrected 

by available experimental energy levels. 



Line fluxes (in unit of xlO 4 phot.cm 2 s x ) 



Index 


Aobs(A) 


Athco(A-) 


Ion 




Transition 




Procyon 


a Cen B 


a Cen A 


e Eri 


1 


43.743 


43.757 


Si XI 


2s 2 




2s3p 


'Pi 


0.54±0.08 


0.35±0.06 


0.23±0.05 


0.33±0.05 


2a 


46.283 


46.294 


Si XI 


2s2p 


3 P1- 


2s3d 


3 D 2 


0.25±0.07 


0.19±0.05 


0.24±0.06 


0.15±0.04 


2b 




46.260 


Si XI 


2s2p 


3 Po - 


2s3d 


3 Di 










2c 




46.310 


Si XI 


2s2p 


3 P1- 


2s3d 


3 Di 










3a 


46.391 


46.401 


Si XI 


2s2p 


3 P 2 - 


2s3d 


*D 3 


0.40±0.08 


0.19±0.05 


0.39±0.07 


0.24±0.05 


3b 




46.404 


Si XI 


2s2p 


3 P 2 - 


2s3d 


3 D 2 










4 


48.720 


48.730 


Ar IX 


2p 5 3s 'Pi 


- 2p 6 


'So 


0.23±0.06 




0.14±0.07 


0.13±0.04 


5 


49.207 


49.229 


Si XI 


2s2p 


!P1- 


2s3d 


l D 2 


1.44±0.14 


0.67±0.07 


0.64±0.07 


0.71±0.06 


6 


52.306 


52.295 


Si XI 


2s2p 


'Pl- 


- 2s3s 


'So 


0.75±0.11 


0.43±0.06 


0.39±0.07 


0.40±0.05 



Capella 



0.61±0.07 
0.60±0.08 



0.43±0.07 

0.16±0.05 
1.33±0.09 
0.73±0.08 



Notes: The same label with different lowercase indices (e.g., 3a and 3b) indicates blended line. 



Table 2: Diagnosed electron temperature in Si XI emitting region for 5 stars. 



Stars 



Ri 



T e (MK) 



R2 



% (MK) 



Procyon 
a Cen B 
a Cen A 
e Eri 
Capella 



1.39 ±0.41 
1.24 ±0.32 
1.72 ±0.57 



1.65 
2.10 
< 2.47 



1.79 

0.80 
1.80 
0.88 



1.20 ± 0.33 2.30;?; 10 



01 



1.20 ± 0.27 2.30^ 



0.64 ±0.14 
0.72 ±0.20 
0.76 ±0.24 
0.62 ±0.13 
0.67 ±0.12 



2.24^ 



+5.55 

0.85 
+4.17 
0.69 



1.52 

1.23 

2.76l L il 
2-02±g : £ 



Table 3: Diagnosed electron density (in unit of xlO 9 cm 3 ) in Si XI emitting region for our 
sample. 



Stars 


#3 


n e 


i?4 


n e 


Procyon 


2.95 ±0.86 


< 87.86 


1.88 ±0.65 


> 0.03 


a Cen B 


3.16 ±0.88 


< 5.00 


2.26 ±0.32 


< 0.95 


a Cen A 


2.13 ±0.82 


> 0.52 


1.62 ±0.52 


> 0.44 


e Eri 


2.46 ±0.73 


> 0.25 


1.65 ±0.55 


> 0.42 


Capella 


2.78 ±0.67 


> 0.06 


1.71 ±0.51 


> 0.33 



-16- 



Table 4: Present results of temperature and those derived from H- and He-like N ions (Ness 
et al. 2002). Lya/He represents the flux ratio between the Lja line and the sum of triplet 
lines of N VI. 







T e (MK) 




Stars 


Present 


Lya/r Lya/He 


G = i±L 

r 



Procyon 


1 65 +L79 

l - uo -0.80 


1.09 ±0.13 


1.41 ±0.07 


1.28 ±0.46 


a Cen B 


2 10 +L80 

z - iu -0.88 


1.79 ±0.16 


1.43 ±0.22 


< 1.8 


a Cen A 


< 2.47 


1.35 ±0.09 


1.15 ±0.09 


< 1.2 


e Eri 


2.301™ 


2.47 ±0.18 


1.95 ±0.30 


< 2.4 


Capella 


9 qn +1 - 57 

Z - OU -0.87 


2.73 ±0.08 


2.02 ±0.10 


0.46 ±0.28 



Table 5: Present results of electron density and those derived from the ratio between % and 
/ lines of He-like N VI ion (Ness et al. 2002). 





log 


n c /cm 3 


Stars 


Present 




Procyon 


> 7.48 


9.96 ±0.23 


a Cen B 


< 8.98 


9.99 ±0.65 


a Cen A 


> 8.64 


9.95 ±0.30 


e Eri 


> 8.62 


10.35 ±0.33 


Capella 


> 8.52 


9.86 ±0.12 



Table 6: Filling factor and electron pressures p derived from n e corresponding to the best-fit 
value of R4. For comparison, the filling factor from O VII are also listed (Ness et al. 2004). 





P 


filling factor / 


in % 


Stars 


dyn/cm 2 


Pres. 


APEC 


O VII 


Procyon 


2.69 


3.78 


2.84 


1.31 


a Cen B 


0.70 


1.44 


1.09 


0.62 


a Cen A 


5.71 


0.14 


0.10 


0.18 


e Eri 


10.86 


0.25 


0.19 


0.62 


Capella 


5.57 


0.45 


0.34 


0.64 
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Table 7: Real constraints on the pressure p (in dyn/cm 2 ) and the filling factor / (in %). 





Procyon 


a Cen B 


a Cen A 


e Eri 


Capella 


p 


> 0.07 


< 2.75 


> 0.57 


> 1.33 


> 1.05 


f 




> 0.36 


< 1.40 


< 2.04 


< 3.30 
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(a) Energy level 
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(b) Weighted oscillator strength gf 
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Fig. 1. — Comparison of our calculated atomic data with other available data, (a) Compar- 
ison of different theoretical energy levels with available experimental ones (in unit of Ryd). 
x-axis is the experimental energy from NIST (open symbols) or the results of Coutinho & 
Trigueiros (2001) (filled symbols), while y-axis represents ratios of different calculation vs 
experimental one. Different symbols correspond to the ratio of results of different calcula- 
tions vs experimental ones, o: present results, o: theoretical data in the CHIANTI database. 
The solid horizontal lines represents 2% uncertainty range, (b, c) Comparison of our results 
with the input data of CHIANTI code, x-axis is the present result, while y-axis denotes the 
ratio of the input data of the CHIANTI code vs present one. Filled up-triangle denotes the 
transitions interested in this work. 
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1.0 L= — — „ — ' — ' L 

10 10 10 10 10 10 

Electron density n e (cm 3 ) 

Fig. 2. — The theoretical line ratios of Si XL Upper, plotted as a function of the electron 
temperature at n c =1.0xl0 9 cm~ 3 . Lower, plotted as a function of the electron density, solid 
line corresponding to logT c (K)=6.2. 
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Electron Density n e (cm 3 ) 

Fig. 3. — The theoretical Si XI emission line ratios vs the electron density. The solid line 
corresponding to logT c (K)=6.2. 



0.03 




43.4 43.6 43.8 44 48.8 49 49.2 49.4 

Wavelength (A) Wavelength (A) 



Fig. 4. — Spectrum of Procyon observed with LETGS and best-fit spectra (smooth lines) 
in wavelength ranges of 43.4 — 44.0 A and 48.7 — 49.5 A. Dot lines represent the continuum 
levels. Prominent lines have been labelled above curves. 
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5.5 6.0 6.5 7.0 

logT (K) 

Fig. 5. — Comparison between the theoretical line intensity ratios and R2) and observed 
ratios of five stars. Opened symbols with error bars are observed values of Ri, while filled 
symbols are values of i?2- 



3.7 




2 1 ' ' 1 1 1.1.1 1 1 

' 10 7 10" 10' 10 1 " 10" 



Electron Density (cm 3 ) 



Fig. 6. — Comparison between the theoretical line intensity ratios (R3 and R4) and observed 
ratios of our sample. Opened symbols with error bars represent observed values of R3, while 
filled symbols are for R4. 



